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Abstract 

We compare the predictions of a range of existing models based on the Vector 
Meson Dominance hypothesis with data on e~^e~ vr+vr" and e~^e~ f-^f-^ 
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cross-sections and the phase and near-threshold behavior of the timehke pion form 
factor, with the aim of determining which (if any) of these models is capable of 
providing an accurate representation of the full range of experimental data. We 
find that, of the models considered, only that proposed by Bando et al. is able 
to consistently account for all information, provided one allows its parameter a to 
vary from the usual value of 2 to 2.4. Our fit with this model gives a point-like 
coupling 77r'^7r~ of magnitude 2± — e/6, while the common formulation of VMD 
excludes such a term. The resulting values for the p mass and tt+vt" and e~^e~ 
partial widths as well as the branching ratio for the decay to 'k~^tt~ obtained 
within the context of this model are consistent with previous results. 
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1 Introduction. 



Our aim is to study the various ways to describe the p meson in order to find an 
optimum modelhng able to account most precisely for the known features of the physics 
involving this meson. This is mainly motivated by the fact that a precise knowledge 
of its properties is of fundamental importance in several fields of particle physics. It is 
important to emphasise that our philosophy is to look for the simplest models as these 
are the most useful in application to other systems, due to their ease of implementation. 
Naturally, such models should, as much as possible, respect basic general principles such 
as gauge invariance and unitarity. One must keep in mind that any parameters quoted for 
a given model are relevant only to that model. Indeed, a study of the model-dependence 
of resonance parameters is one of the principal goals of this work. It should also be noted 
that while each of these models is related to some underlying effective field theory through 
an effective Lagrangian, the models we are using here are simple amplitudes arising from 
an assumption of almost complete s-channel resonance saturation. The appropriateness 
of this assumption away from the resonance region can only be judged by quantitative 
studies of higher order (e.g., loop) effects in the corresponding effective field theories. 
While this is clearly a very important task, it is not our concern here and will not be 
considered further. 

For this purpose, we study the strong interaction corrections to one-photon mediated 
processes in the low energy region where QCD is non-perturbative. To do this we shall 
look at two related processes, e+e~ — ^ tt+tt^ and e+e" — > . The effect of the 

strong interaction is obvious in the first reaction and provides a large enhancement to 
the production of pions in the vector meson resonance region [||-^. This enhancement, 
relative to what would be expected for a structureless, pointlike pion, is refiected in 
the deviation of the pion form factor, F.,r{q'^), from 1, and is primarily associated with 
the p meson (where is the four momentum of the virtual photon). This form factor is 
successfully modelled in the intermediate energy region using the vector meson dominance 
(VMD) model 0. VMD assumes that the photon interacts with physical hadrons through 
vector mesons and it is these mesons that give rise to the enhancement, through their 
resonant (possessing a complex pole) propagators of the form 



where my and Ty are the (real valued) mass and the momentum-dependent width. (Here 
we have included only that part of the propagator which survives when coupled to con- 
served currents.) 

Traditionally, VMD assumes that all photon-hadron coupling is mediated by vector 
mesons. However, from an empirical point of view, one has the freedom, motivated by 
Chiral Perturbation Theory (ChPT) to include other contributions to such interactions. 
For instance, in a fit to F-j^lq"^) 0, it was shown that a non- resonant photon-hadron 
coupling can be accommodated merely by shifting the mass and width of the p° by about 
10 MeV. Thus, the values extracted from -F^(g^) for the p mass, nip, and width, Tp, 
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are model-dependent and in quoting values for them, the model used should be clearly 
stated. We note in passing that to define vector meson masses and widths in a process- 
independent way, one should refer to the location of the corresponding complex pole in 
the S-matrix. One should, however, bear in mind that alternate, and more traditional 
definitions of the mass and width not tied to the location of the S-matrix pole (for 
example a definition of as that value of for which the P-wave vrvr phase shift 
passes through 90°) are in general specific to the process employed in the definition. The 
process-dependence of such alternate mass and width definitions has in fact led a number 
of authors to advocate using the S-matrix pole position to provide a process-independent 
definition of the Z° mass and width in the Standard Model 0. 

Naturally, VMD can be applied to many other systems. We can consider the process 
7]' /t] —* 7r+7r~7 taking place through a combination of resonant (such as t]'/-)] ^ p7 — >■ 
TT+TT^^) and non-resonant channels. In this manner an acceptable fit to data can be 
achieved with a range of combinations for the p parameters and non-resonant terms P,^] . 
Recent interest in this process has centred on the non-resonant term, which, if it arises 
from anomalous box and triangle diagrams, provides a possible test of QCD P-[TT 



However, to determine the size of any non-resonant contribution, the resonant meson 
parameters need to be well fixed and thus F-^^q"^) well understood. This is one of 

the main aims of this paper. 

We thus turn our attention to e~^e~ — fi'^fi^ . In modelling the strong interaction 
correction to the photon propagator, VMD assumes that the strong interaction con- 
tribution is saturated by the spectrum of vector meson resonances [0. Therefore, in 
principle, we can extract information on the vector meson parameters (independently of 
the e^e~ tt+tt^ fit) without having to worry about non-resonant processes. However, 
as the vector mesons enter in e~^e~ fi'^fi" with an extra factor of a compared with 
e~^e~ TT+TT" , their contributions are considerably suppressed, making their extraction 
difficult. For this reason, we shall perform a simultaneous fit to both sets of data, in order 
to impose the best possible constraint on the vector meson parameters, and see if existing 
muon data are already precise enough in order to constrain the parametrisation. 

Another way to constrain the descriptions of the p° meson is to compare the strong 
interaction vrvr phase obtained using the various VMD parametrisations determined in 



fitting e"^e~ tt+tt" with the corresponding phase |T3[ obtained using tttt scattering 
data and the general principles of quantum field theory, as well as the near-threshold 
predictions of ChPT. This happens to be more fruitful and conclusive in showing how 
VMD should be dealt with in order to reach an agreement with a large set of data and 
with the basic principles of quantum field theory. 

The hadronic dressing of the photon propagator (the one-photon irreducible self- 
energy H^'^'^(g^)) is also of interest for another reason. The anomalous magnetic moment 



of the muon can now be measured to such accuracy |]I4| that the strong interaction 



correction is important |]I5[. This needs to be completely understood if one is to look for 
physics beyond the Standard Model in this quantity. At present the correction is inferred 
from a{e~^e~ — >■ hadrons) using dispersion theory and the optical theorem or estimated 
with hadronic models (which, being non-perturbative, are difficult to use). The process 
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e"'"e~ — >■ jji'^jj,^ allows for a direct examination of the strong interaction modification to 
the photon propagator. Ideally we would obtain the low energy corrections to the photon 
propagator experimentally, and be able to use QCD perturbatively for higher energies 
(though the threshold above which we can ignore non-perturbative effects is difficult to 
determine |jl6[)- The meson has already been seen in e^e~ — fi~^fi~ IH,^- It is only 
noticeable around the pole region where, due to its small width (4.4 MeV), it produces a 
sharp peak easily seen in the available data. The large width of the p further suppresses 
the p and uj contributions. 

The outline of the present work is as follows. In section |^ we describe the various 
formulations of the VMD assumption and the unitarisation procedure; we also discuss the 
phase definition relevant for our purpose. The fit procedure is sketched in section ^ and 
implemented in section Comparison with the isospin 1 P-wave vrvr phase shift deduced 
from vrvr scattering theory is the subject of section ^, while in section ^ we give the near- 
threshold parameter values that are deduced from our fits of the e~^e~ — > tt+tt^ cross 
section and they are compared to predictions and to other experimental determinations. 
The results obtained are discussed in section |^ where we present our optimal fit for the 
p° parameters and the model which accounts for its properties the most appropriately. 
Finally, we summarise our conclusions in section p. 



2 Vector meson models 



We shall now provide a description of the various models we will use to fit the data for 
both e~^e~ — >■ 7r"'"7r~ and e~^e~ — >■ p^ p~ . The cross-section for e^e~ tt+tt" is given by 
(neglecting the electron mass) 



a 



(2) 



where the form factor, Fj^lq"^) is determined by the specific model. Similarly, F^j_{q^) is 



defined to be the form factor for the muon, and the cross-section for e^e 
given by. 



p^p IS 



a 



3g2 \ 




(3) 



It is worth noting that these standard definitions of F-j^lq"^) and F^{q'^) contain all non- 
perturbative effects, including for example the photon vacuum polarisation, since Eqs. (0) 
and d^) are written assuming a perturbative photon propagator. 

We shall use VMD, the Hidden Local Symmetry model of Ref. [0 (hereafter ref- 
ered to as HLS) and what we will refer to as the WCCZW model (a phenomenological 
modification of the general framework of Ref. [|18|, based on work by Birse as 
well as modifications of these models, which cover or underlie a large class of effective 
Lagrangians describing the interactions of photons, leptons and pseudoscalar and vector 



mesons. Birse [19| has shown that typical effective theories involving vector mesons based 
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on a Lagrangian approach, such as "massive Yang-Mills" and "hidden gauge" (i.e., HLS- 
type) are equivalent. We therefore consider our following examination to be reasonably 
comprehensive. Numerical approaches to the strong interaction in the vector meson en- 
ergy region , which are not based on an effective Lagrangian involving meson degrees 
of freedom do not yet possess the required calculational accuracy for our task. 



2.1 VMD 

The simplest model is VMD itself. As has been discussed in detail elsewhere P,^^ 
VMD has two equivalent formulations, which we shall call VMDl and VMD2. The VMDl 
model has a momentum-dependent coupling between the photon and the vector mesons 
and a direct coupling of the photon to the hadronic final state. The resulting form factor 
is (to leading order in isospin violation and a = e^/Air): 

pVMDl/ 2\ ^ _ VMDl/ 2\ dpirn 

' ^[g2-m2 + ^m,r,(g2)] 

-slT (^') n 2^ V(2^^ ^^''' ■ (4) 

The uj enters into the isospin 1 e^e'^ tt^tt^ interaction with an attenuation factor 



specified by the pure real A and the Orsay phase, ||2^; these can be extracted from 
experiment^ For VMD2 we have. 



pVMD2/ 2\ ^ _ VMD2 dp-KTv VMD2 ^ a i<t,2 (^\ 



The form factor for the muon, however, in both representations (i=l, 2) is given by 

^VMDi ^ 1 + ^e^[^VMDi(^2)p^ \ I (g) 

V q^ - mf, + inivTviq^) q^ 

which is consistent with previous expressions for the 0-meson In higher order {i.e., 

in all but the minimal VMD picture) there will also be contributions from non-resonant 
processes (such as two-pion loops), but these are expected to be small near resonance and 
the non-resonant background is, in any case, fitted in extractions of resonance parameters 
from the experimental data. The photon-meson coupling, egy-y is fixed in VMD2 by 

Ty^e^e- - ^9vr (7) 
The (dimensionless) universality coupling, gy, is then defined by 



V 



■^Note in Ref. |Q that A and (f) were defined through the S matrix pole positions (equivalent to Eg. (p|) 
with constant widths). In our fit procedure, A is connected with the width r(w — > 7r"'"7r~) (see Ref. M). 
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for VMD2 PT],pB[. This coupling (and universality) has been most closely studied for 
the p meson. A gauge-like argument P, ^ suggests that the p couples to all hadrons 
with the same strength gy (universality) However, experimentally, universality is 



observed to be not quite exact |]26|, so we introduce the quantity e (to be fitted) through 

+ e) (9) 



VMD2 
Up-y 



9p-n 



where gp-,^-,^ and e are to be extracted from the fit to e tt+tt^ . For VMDl, it can 

be seen that the photon-meson coupling results from replacing the mass term in Eq. (P) 
by [iill^- In this case Eq (^ should be replaced by 

^'-(1 + e) . (10) 



^VMDl 
9V"i 



Qp-n 



One can easily see that in VMDl the hadronic correction to the photon propagator 
goes like and so maintains the photon pole at = 0. For VMD2, it is not obvious 
that gauge invariance is maintained until one considers the inclusion of a bare photon 
mass term in the VMD Lagrangian that exactly cancels the hadronic correction at = 0. 
This argument assumes a p — 7 coupling of the form em? / gp and a bare photon mass 



{e'^m'^p/ gp) (the calculation is presented in detail in Ref. f^). In the presence of a finite 
e, as in Eq. (0), gauge invariance is similarly preserved by including a photon mass term 



1 + / g'l) leading to a massless photon as expected 



(e^m2( 

The presence of a finite e does affect the charge normalisation condition ^,^(0) = 1, 
but this is merely an artifact of the simple p propagators we are considering. A more 
sophisticated version, such as used by Gounaris and Sakurai which fully accounts 
for below threshold behaviour, maintains Ft^{Q) = 1 in the presence of e 7^ 0. One 
could, alternatively, include an s dependence to e such that e(0) = 0. In any case, the 
phenomenological significance for the physical p (for the data region we are fitting) is 
negligible, we achieve an excellent fit to the data with the simple form we use but are 
not advocating its use outside this fitting region, namely above the two-pion threshold. 

The choice of the detailed form of the momentum-dependent width, Ty{q^), allows one 
certain amount of freedom. As the complex poles of the amplitude are field-choice and 
process- independent properties of the S-matrix [31| one could, on the one hand, expand 
the propagator as a Laurent series in which the non-pole terms go into the background 
3^ . Alternatively, one could use the Z— wave momentum-dependent width to account 
for the branch point structure of the propagator above threshold (g^ = 4m^) P,pU|,P^. 
This form for the momentum dependent width arises naturally from the dressing of the 
p propagator in an appropriate Lagrangian based model (see Sec. 5.1 of Klingl et al. for 
a detailed treatment EH]) and, for such models, is given by 



Tp{q' 



A/2 



'11^ 



introducing the fitting parameters T p (the width of the p° meson at q^ 
which generalises the usual /-wave expression H] to model the fall-off of the p mass 



rrip) and A, 
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distribution; the usual case {i.e. A = 1) is associated with a p couphng to pions of the 
form 5'p7r7rP^(vr+9^7r~ — 7r~(9^7r+), with Qpt^t, independent of g^, as shown in Ref. p6[. Note 



that the parameters A, nip and T p are model-dependent (as we see in the tables of results, 
different models yield different values for these parameters). Note that in Eq. ( [Tl| ) we 
have defined the pion momentum in the centre of mass system 

pM^) = yq^-Aml (12) 

Before closing this section, let us remark that the pion form factor associated with 
VMDl (Eq. (ID) fulfills automatically the condition -F^(O) = 1 whatever the value of the 
universality violating parameter e (see Eq. (PD])). This is not the case for the pion form 
factor associated with VMD2 (Eq. (H)), as can be seen from Eqs. (^ and (P). Here, 
we will concern ourselves exclusively with fitting data in the above threshold region 
and simply note it is a relatively straightforward matter to generalise the VMD models 
considered to satisfy this condition. Detailed considerations of this issue are left for future 
work. 

While of course VMDl and VMD2 are equivalent in the limit of exact universality 
if one keeps all diagrams, (i.e., if one works to infinite order in perturbation theory), in 
any practical calculation one cannot do that and so, in practice, these two expressions 
of VMD can give different predictions in general, even if exact universality is imposed. 
Moreover, if one releases (as we do) this last constraint, equivalence of VMDl and VMD2 
is not guaranted, even in principle. 



2.2 The HLS Model 



The Hidden Local Symmetry (HLS) model |^,|T3 introduces a parameter a for the p meson 
within a dynamical symmetry breaking model framework. This a relates the constant 
Qp-^T, to the universality coupling Qp via 



QpTTTT 



2 



The resulting form factor for the pion is 



9py 



QpTTTT 



1 



-Ae 



(13) 



(14) 



2 ' (q'^ - ml + imJ'niQ'^)) ' q"^ - ml + im^T^{q'^) 

The original HLS model preserved isospin symmetry and so did not include the u. Isospin 
breaking has recently been studied in a generalisation of the HLS model 0], however, 
here we have for simplicity employed the same u terms as used for VMD. The relations 
equivalent to Eqs. (S) and (pi) for the p meson are now 



a 

2 gp-iv: 



m. 



9p 



(15) 
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We see that setting a = 2 reproduces VMD2 in the hmit of exact universahty. However, 
we wish to keep a as a free parameter which we can fit to the data. Note that in the 
HLS model universahty violation and the existence of a non-resonant coupling 77r+7r^ 
are related. Note also that universality violation can be introduced in the HLS model 
without violating the constraint -F7r(0) = 1, in a natural way. 

The muon form factor for the HLS model is exactly the same as for VMD2 (see 
Eqs. d) and (|15D). 



2.3 WCCWZ Lagrangian 



Birse has recently discussed the pion form factor arising from the WCCWZ La- 
grangian in which the vector and axial vector fields transform homogeneously under 
non-linear chiral symmetry. The scheme imposes no constraints on the couplings of the 
spin 1 particles beyond those of approximate chiral symmetry. Birse's version of the form 
factor is (isospin violation is not considered) 



9ifi r , f'. 



fi q2 _^2^ impTp{q^) 

where the first two terms on the RHS are those arising from the WCCWZ Lagrangian. 
The q'^ piece grows at large q^ in a way incompatible with QCD predictions (for a discus- 
sion of matching the asymptotic prediction to a low energy model see Geshkenbein p5|). 
The /2 contribution has been added by Birse to modify this high energy behaviour toward 
that expected in QCD. To this end, Birse sets 

/2 = gifi = 4 (17) 

nip 

and recovers the universality limit of the form factor in which VMDl and VMD2 are 
equivalent (in the zero width approximation). Note that a g^-dependence of the non- 
resonant background is what one would, in general, obtain from the WCCWZ framework, 
implemented in its most general form, which relies only on the symmetries of QCD. In 
constructing a phenomenological implementation we have, however, simplified the most 
general form, adding what amounts to a minimal g^-dependence to the background term 
of VMDl. The resulting form factor, which we will refer to as the WCCWZ model, is 
then 



WCCWZ/ 2X . i0 

F^(q^) = l + bq^ - - -— 



where we keep b as an independent parameter to be fit. 

The WCCWZ model, thus, has one more free parameter than VMDl. We have added 
the uj contribution as above for all other models. The muon form factor is exactly the 
same as for VMDl (see Eq. @). 

It is important to note that the WCCWZ model allows one to have a non-resonant 
term which can be mass dependent, cf., VMD2 which carries only resonant contributions 
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or VMDl or the HLS models which both exhibit only constant non-resonant contributions 
to the pion form factor. The expression in Eq. (|18]) exhibits an unphysical high energy 
behaviour; however, as we are only interested in the pion form factor at low energies 
(below 1 GeV), this feature is not relevant. Of course, one can consider that such a 
polynomial structure at low energies represents an approximation in the resonance region 
to a function going to zero at high energies 0. 

The model VMD2 contains only resonant contributions, whereas VMDl, HLS and 
WCCWZ also contain a non-resonant part. For VMDl and HLS, this term is constant 
{i.e., pointlike) as in standard lowest order QED. We shall frequently refer to this term as 
a direct 'yirn contribution or coupling. In the case of WCCWZ, this non-resonant term 
also contains a g^-dependent piece which clearly indicates a departure from a point-like 
coupling; nevertheless, we shall also refer to it as a direct jtitt coupling for convenience. 



2.4 Elastic Unitarity 

When A 7^ 1 the pion form factor described above in the VMD, HLS and WCCWZ models 
does not in general exactly fulfill unitarity. This is because one employs simultaneously 
a bare, undressed p^n^n^ coupling (given by the tree-level couphng, gpTrn) and a dressed 
p° propagator, as signalled by the presence of a momentum-dependent width in Eq. ([Tl|) 



(see, e.g., Ref. @]). To see why this generally creates a problem with unitarity, consider 
the TTTT —>■ TTTT P-wave 1 = 1 strong interaction amplitude. It is most convenient here to 
employ the N/D formalism for the amplitude T (see, e.g., Ref. \^% ). 



It has been known for some time that the J = / = 1 vrvr amplitude is consistent with 
being purely elastic from threshold up to ~ 1 GeV [p^ -|40[]. Therefore, in this region, we 
can write0 

S = 1 + 2ip^T = exp [2iSl] (19) 

from which it follows that 

where has been defined in Eq. (0). In the region where the p° meson essentially 
saturates the J = / = 1 vrvr wave, the D function may be approximated by the inverse p° 
propagator (s = g^), 

Dis) = s-ml + tmpTp{s) . (21) 

D in Eq. (pp]), is an analytic function of s, having as sole singularity in the physical 
sheet a cut along the real axis (s > 4m^). Correspondingly, the singularities of the N 
function on the physical sheet are all located on the real axis at negative values of s (for 
example, singularities produced by exchanges in the t and u channels projected out onto 

''For instance, 1 + bq^ can be considered as the first terms of the Taylor expansion of a function like 
1/(1 — bq^) as suggested by [|6). 
^The notation here is Sf . 
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the P-wave). Moreover is real on the real axis above threshold (4m^). It should also 
be noted that the and D functions in Rel. (0), are defined up to a multiplicative 
arbitrary function /(s), meromorphic in the complex s-plane. The choice used in Eq. (|2l| ) 
corresponds to a particular choice of /(s). Put simply, the content of Eq. ( [21| ) is that the 
phase of D is well approximated by the phase of the propagator in the vector meson 
resonance region. 

As a consequence of unitarity and analyticity, is connected with the discontinuity 
of D across the physical region (s > 4m^) through 

lim \D{s-ie)- D{s + ie)]=2ip^N{s), (22) 

which gives 0, 

N{s) = (23) 

As the N function is closely connected with the pt:^t:~ vertex, this last relation implies 
some dressing of the vertex coupling. This relation illustrates the effect of the parameter 
A, whose role is simply to model the contributions from the left-hand singularities of 
the scattering amplitude. Then using Eqs. (^ - ^3|), it is easy to check that \S\ = 1 is 
automatically satisfied in the physical region below ~ 1 GeV, as required by unitarity, 
and that 

tan 51 = !^^. (24) 
— s 

Thus from Eq. (^) we can conclude that 6l can be well approximated by the negative 
of the p propagator phase, if one considers resonant contributions only. 

One can then show |0 that unitarity can be restored to the model treatments above 
by replacing the coupling gpj,.,, in Eqs. (^, (|^), ([l^) and ( [TSl) with 



GpiQ') = ^|Q^^^^P^.^n-iq') (25) 



where q = y/q^. This replacement leads to the unitarised versions of our VMD models. 
The connection between this "dressed" vertex function and gpmr gives (cf., Eq. (4.16) of 

Ref. m 



QpTTTV 

from which we see that 



6^^7^r„ (26) 



fl-pTTTT = Gp{ml). (27) 

It should be noted that the left hand side of Eq. (p5[ ) becomes constant - and then 
coincides with Qpt^t^ - if and only if A = 1 (see Eq. (|ll|)). Therefore, if A = 1 Eqs. §), (D, 
(|I^) and (|TS|) are already unitarised and the singularities of (see Eq. (P^) are simply 
a branch point going from s = to — oo, where s = is also a pole. 
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Strictly speaking, one should similarly unitarise the uj contribution. However, since 
the uj is very narrow, its contributions are significant only over a very limited range of 
. Hence, it is sufficient to employ a constant cuvrvr coupling (or equivalently, constant 
A) in Eqs. (^, (^), (|T^) and (p^S]). The value of A may be determined from the branching 



fraction of — vr+vr^, as was done in Ref. or fitted. Therefore, the only actual free 
parameter in the lo contribution is the Orsay phase, 0. This will in no way affect the 
description of the p° itself, which is our main concern. Finally, unitarisation clearly does 
not affect the expressions for the muon form factor. 

The unitarised version of VMD2 coincides with the phenomenological model called 
M2 in Refs. Here, however, we shall fit gp^ and then BR(p° e"'"e~), which was 

previously simply fixed at its PDG |^T| value 0. The unitarised version of the HLS model 



is close to the phenomenological model M3 of Ref. 

We note that the non-resonant term introduced by VMDl and HLS models is constant 
at leading order whereas Mi of Refs. P,^] has a g^-dependent non-resonant piece. In 
this respect the WCCWZ model is close to Mi. 

2.5 Phase of FT^{q^) and Phase of the tttt tttt AmpUtude 

From general properties of field theory (mainly, unitarity and T-invariance) , it can 



be shown [42] that for s = q real above threshold, we have 



F^{s)=exp[2z6l] F;(s) , (28) 
up to the first open inelastic threshold. Using Eq. (|19|) and Eq. ( |20| ) - which defines the 



N/ D formalism-, this relation gives 



D( 



(29) 



where D{s) has been defined by its general properties in the preceding section. In order 
that F^{0) = 1, the D function should fulfillQ -D(O) = 1. The functions D appearing here 
and in Eq (pOD can be chosen identical without any loss of generality. It is obvious from 
Eq. (p9|) that -F^(s) carries the same phase as 1/D{s) and as the tttt tttt amplitude. 

The identification of the function D in Eq. ( pO]) with a resonance propagator, namely 
the p° meson (see Eq. (^I])), would be motivated in the present case| by the fact that the 
amplitude is dominated by a single resonance (the p° meson itself) P^^D|. This does 
not mean that higher mass resonances, which surely exist 0,0, have a (strictly) zero 
contribution in our mass range, but simply that their magnitude is negligible compared 

^This parameter was taken as fixed in Refs. |^,|^ because tlie autliors believed tliat the PDG value for 
Tili{p° — *■ e+e") relied on e'^e^ data and other data. As there is no other data than e+e^ annihilations 
for this branching fraction, it should be fitted. 

^ In phenomenological applications, one could prefer requiring a condition on a physically accessible 
invariant mass region as, for instance, the two-pion threshold, where we have i^7r(4?TiJ) — 1.17 ± 0.01 
from ChPT estimates, as it will be seen below. 

^Namely: s < 1 GeV^, I = I = 1 tttt tttt partial wave. 
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to that of the p meson. In this case, the Breit-Wigner parametrisation is flexible enough 
in order to absorb these small contributions into a correspondingly small change of the 
parameter values with respect to their (unknown) "true" values. This is also valid for a 
possible small non- resonant hadronic contribution to the amplitude. 

As far as hadronic contributions to the pion form factor are concerned, there are also 
contributions from the uo and mesons. Because of their small width, we can add their 
contributions to Ft^{s)\ this will produce departures from Eq. (^), however always in 
very limited invariant mass interval. Moreover, they do not contribute to the 1 = 1 = 1 
7171 phase shift for obvious reasons. 

Putting aside the question of u and mesons for the reasons just given, it remains 
to recall that there is a direct 'jtttt contribution to the form factor present in the models 
VMDl, HLS, and WCCWZ, whereas VMD2 has no such contribution. These additional 
direct contributions, which are constant for VMDl and HLS models, modify the total 
phase near threshold where the hadronic (p^) contribution is small in magnitude. There- 
fore, it is appropriate to take the full phase of F^{s), rather than the phase of the p° 
propagator only, which allows us to extract the exact behaviour of the vrvr — tttt phase 
{6l) in the threshold region. 



3 Simultaneous Fits of e+e n^n and Data 

Fitting the e~^e~ data from threshold to about 1 GeV involves three well-known res- 
onances, p°, uj and (f). As the last two are narrow, their parametrisation is relatively 
simple. But, due to its broadness, the p° meson has given rise to long standing problems 
of parametrisation (see Refs. p,^^] and previous references quoted therein). Moreover, 
as we have mentioned, one can ask whether experimental data require the existence of a 
non-resonant 77r+7r^ coupling. The conclusion of Refs. is that data on F-j^lq"^) alone 
are insufficient to answer this question. This is also relevant to the test of QCD proposed 
by Chanowitz [|,|ig. 

In addition to the mass and width, in the context of the class of models having widths 
of the form given in Eq. (pJ]), one requires 



only one additional parameter. A, to define the p shape. The resulting fit turns out 
to depend not only on the non-resonant coupling, but also on whether unitarisation is 
used or not. One approach |Q to this problem is to perform a simultaneous fit of all 



e~^e~ — 7r~^7T~ data and e~^e~ — > p~^p~ data [^. Indeed, if the data are precise 
enough, we could see a non-resonant coupling in e'^e^ tt+tt^, which will (of course) 
be small in e~^e^ fi^p~. Therefore, from first principles, a simultaneous fit to both 
data allows us to decouple the p from any non-resonant 7'7r^7r~ coupling. Naturally, 
to be of any use in this, the e~^e~ p~^p~ data would have to be very good. Until 
recently relatively precise measurements were available only for the region around the 



mass []l],§]. However, a new data set collected by the olya collaboration is available 
and covers a large invariant mass interval from 0.65 GeV up to 1.4 GeV. We shall see 
shortly whether it is precise enough to constrain the p parameters. Thus, the data sets 
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which will be used for our fits are those collected by DMl, olya and CMD which are 



tabulated in (for e~^e vr+vr ) and only the olya data of |^ for e~^e — >• /i"*"// ; 
these data sets do not cover the peak region. 

4 Results of e^e~ data analysis 

In all of the previously described models, except for WCCWZ, the fit to e+e" — > vr+vr" 
and e^e~ — *• ix^ ii~ data depends on only five parameters. The first four are the three p 
meson parameters {mp, T p and A) and the Orsay phase (0). These are common to both 
VMD and HLS. In the VMD models an additional parameter, e, has been introduced 
in order to account for universality violation (see Eqs. (§) and (|ToD), while in the HLS 
model this parameter is replaced by a (see Eq. (0)). These last parameters allow us to fit 
the branching fraction p° — > e^e~ within each model in a consistent way. The WCCWZ 
model depends on one additional parameter, 6, which permits a more flexible form for 
the non-resonant contribution, as compared with the VMDl or HLS models. Let us note 
that introducing e^e~ p^p~ in our fit procedure together with e~^e^ tt+tt^ does 
not require further free parameters. 

Generally speaking, the parameter named A in the VMD models above determines 
the branching ratio Br(ct; —>■ 7r"'"7r~) and should be set free since its value is strongly 
influenced by the data on e^e~ tt^tt" we are fitting. However, in order to minimise 
the number of fit parameters at the stage when different models are still considered, we 



fix its value from the corresponding world average value |^T| of Br(ci; tt+tt^). We shall 
set A free for our last fit, in order to get an optimum estimate of Br(co' —>■ tt+tt"); this 
will be done only for the model which survives all selection criteria. 

Finally the fits have been performed for both the standard VMD, HLS and WCCWZ 
models and their unitarised versions, for both e+e~ tt+tt^ data alone and simulta- 
neously with e^e~ — fi~^fi~. As all measurements in the region of the resonance for 
each of these final states are not published as cross sections they are not taken into 
account in our fits. When fitting the data, we take into account the statistical errors 
given in for each e~^e~ tt+tt^ data set. DMl and CMD claim negligible systematic 
errors (2.2% for DMl and 2% for CMD, while the statistical errors are typically 6% or 
greater); these errors can thus be neglected with respect to the quoted statistical errors. 
OLYA claims smaller statistical errors but larger systematic errors: these two errors have 
comparable magnitudes from the p° peak to the mass. We do not expect a dramatic in- 
fiuence from neglecting these systematic errors, except that this would somewhat increase 
the value at minimum and hence worsen slightly the fit quality. 

The results are displayed in Table |I| (non-unitarised models) and in Table || (unitarised 
models). We show the fitted parameters in the upper section of each table, while in the 
lower part we provide the corresponding values for derived parameters of relevance. 
We find that, disappointingly, the new muon data places no practical constraint on 



the p parameters extracted from F^(g^). Indeed, the central values for the fit parameters 
at minimum are practically the same when fitting only tttt than tttt + pp data, the 
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Parameter 



VMDl 



VMD2 



HLS 



PDG 



e 

HLS a 
rup (MeV) 
Tp (MeV) 
(f) (degrees) 
A 

xVdof (tttt) 
xVdof (tttt + /i/i) 



U.ZiU_Q 018 



751.4^18^ 
146.0 ±2.2 



113.8: 



-5.2 



4,49+0-35 



6.9 

.35 

-0.43 

63/77 
105/115 



1 fiS+°-°°'^ 

U.iUO_Q 008 



776.74±2.2 

145.10t?] 

106.3 ±4.5 

1 61+°'3^ 
-■-•O-L-o.ai 

148/77 

194/115 



2.399t°:°l 
755.1111 
143.32t^:g 

r4.6 

-5.7 



769.1 ±0.9 
151.0 ± 2.0 



120.6: 



64/77 
108/115 



9fn (GeV2) 
g., (GeV2) 



r(p 



e ' e 



(keV) 



0.113 ±0.003 



2.91 ±0.05 
6.72l°:i 



0.119 ±0.002 

2.76 ±0.03 
6.74 ±0.20 



0.115 ±0.003 



2.84 ±0.05 
6.341°:^? 



0.120 ± 0.003 
0.036 ± 0.001 

6.77 ±0.32 



Table 1: Results from fits to Fj^^q^) and F^{q^) without unitarisation 
for the VMDl, VMD2, and HLS models. Shown for comparison are the 
Particle Data Group quoted values ET . 



errors become slightly larger in the second case (because of the magnitude of the errors 
in the resonance region with the /i/i final state). The errors quoted in Tables |l] and ^ 
are those obtained when fitting only tttt data up to 1 GeV/c. Thus the single existing 
data set on the muon final state is not precise enough for our purposes. To improve the 
situation we require more accurate muon data below ~ 1 GeV/c. This is illustrated by 
Fig. |I| which shows the muon data together with one of the best fits (namely unitarised 
VMD2). 

Another striking conclusion is that it is generally possible to achieve a very good fit 
to the pion data whichever model is used, unitarised or not; the single exception to this 
being non- unitarised VMD2 (which is the usual model for p physics). Correspondingly, 
the significant model dependence of the extracted p mass should be noted. We do not 
give results in case of the non-unit arised WCCWZ model, as in this case the solution 
converges to 6 ~ and then coincides with VMDl. We do not present pion data curves 
for VMDl and HLS as they show results indistinguishable from Fig. 4 of Ref. which is 
probably the best possible fit (x^/dof=61/77). Fig. |^ shows the fit obtained using non- 
unitarised VMD2; one clearly sees that the model fails to describe the mass region below 
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Parameter 


VMDl 


VMD2 


HLS 


WCCWZ 


e 

HLS a 
WCCWZ h (GeV-2) 
nip (MeV) 
ip (^MeVj 
(degrees) 
A 

y^/dof fTrvr) 
X^/dof (vrvr + /i/i) 


0.167 ±0.008 

774.67 ±0.65 
A Ai A ^ -1-1 fin 

94.7 ±4.3 

i.Ud»_o.o85 

65/77 
104/115 


0.215 ±0.010 
780.37 ±0.65 

IKK /I/I -1-1 OK 

98.8 ±4.4 
0.567 ±0.055 
81/77 
128/115 


2.364 ±0.015 
775.15 ±0.65 

l/IVfiV-l-l /IV 
14 / .0 ( ± I At 

105.1 ±4.3 
1.056 ±0.042 
65/77 
111/115 


0.142 ±0.014 

-0.319l[l:l?? 

770.8911:^ 
1 /in 

101.7 ±5.3 

i.UZO_Q 269 

61/76 
103/114 


^7P7 (GeV2) 
r(p e+e-) (keV) 


0.118 ±0.001 
2.81 ±0.03 
6.70 ±0.11 


0.122 ±0.001 
2.94 ±0.04 
6.99 ±0.16 


0.114 ±0.001 
3.08 ±0.03 
6.23 ±0.11 


0.133 ±0.007 
2.08 ±0.04 
8.62 ±0.46 



Table 2: Results from fits to FT,{q^) and F^^^q^) for the unitarised VMDl, 
VMD2, HLS and WCCWZ models. 



the uj mass0, while the region from the uj to the mass is quite correctly reproduced. 
It is doubtful that taking into account higher mass mesons could cure this problem. On 
the other hand, Fig. ^ shows the fit obtained using unitarised VMD2 which is of good 
quality throughout the mass range. Its quality is better than that of M2 (see Fig. 3 of 
Ref. 0) simply because (actually, e) is treated here as a free parameter, as it should 
be. 

For the non-unitarised HLS model, the p mass and widths appear similar to those 



obtained in the fit by Bernicha et al. P2|, {nip = 757.5 ± 1.5 MeV and 142.5 ± 3.5 



MeV respectively), which used propagators with constant widths defined in terms of 
the complex pole locations. The reason can be traced back to the fact that the HLS 
model mimics quite well a Laurent expansion when there is no mass dependence in the 
numerator of the p contribution as in Eq. ([T^) . The unitarised fits (shown in Table 
however, all show much higher p masses than do the non-unitarised fits. 

The values obtained for the widths r(p —>■ tt+tt^) and r(p —>■ e^e~) and the Orsay 
phase, 0, are in the expected range. The value for gfp^ is approximately three times that 

^The region which strongly worsens the given in Table is the mass interval from 400 to 600 MeV. 
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of g^^ as expected from SU(3) symmetry and ideal mixing and gp^^/An is always at 
the expected value (about 3) except for the unitarised WCCWZ model which provides 
a much smaller value. Moreover, the smallness of the statistical error on gp^ should be 
noted, i.e., the uncertainty is essentially all due to model dependence. 

The value for A deviates from 1 [see the discussion surrounding Eq. (|Tl|) ] quite sub- 
stantially for the standard VMDl and HLS models, although the unitarised versions 
return a more standard value. As A is a priori an effective parameter which accounts for 
higher order effects in the perturbation expansion and/or dynamical left-hand singular- 
ities, this is not yet grounds for rejection of the non-unitarised versions of the models. 
On the other hand, unitarised VMD2 provides, as expected, a value for A close to that 
found using M2 in Ref. 0. We also note that universality (e = or a = 2) is broken 
significantly in all fits at the level of ~ 20%. Finally, the effective parameter b introduced 
in the unitarised WCCWZ model stays relatively small and allows one to obtain a more 
conventional mass for the p° meson |^ . 

It should be noted that, fixing A to 1 in the non-unitarised versions of VMDl and 
HLS, we essentially recover the solutions given in Table 0, as remarked in subsection |2.4| ; 
the value is practically as good. The situation is completely different for VMD2, for 
which A is found significantly different from 1. As there is no general requirement based 
on unitarity and analyticity which constrains the value of A, nothing can be concluded 
from this observation. It is however interesting that unitarisation of the VMDl and HLS 
models returns A ~ 1; comparing for instance the results for the HLS model given in 
Table with the corresponding results for model M3 in Ref. j^, clearly shows that this 
is a consequence of having released T{p^ — ^ e^e^) in the fits. 

As far as one relies only on the statistical quality of fits for the cross section e'^e~ 
7r+7r~, the existing data do not allow one to determine the most suitable way to implement 
vector meson dominance, except to discard the non-unitarised version of VMD2 which is 
clearly disfavoured by the data. On the other hand, the possible values for the mass 
cover a wide mass range: from 750 MeV to 780 MeV. The single firm conclusion which 
can be drawn from the above analysis is that, whichever is the VDM parametrisation 
chosen, unitarised or not, one always observes a small but statistically significant signal 
of universality violation: e ~ 0.20 (instead of 0) or a ~ 2.4 (instead of 2). 

The question, therefore, remains as to whether it is possible to find other criteria 
to distinguish between the various ways of building effective Lagrangians involving the 
vector mesons. 



5 Comparison with tttt Phase shift Analysis 

Until now, we have focussed on the various ways of implementing the VMD assumption 
when describing the e^e~ — > tt+tt" cross section. As is clear from Section 0, the way to 
express VMD is not unique. A representative set of possible models has been given in 
the previous section. We say "representative" in the sense that it is difficult to imagine 
a class of models very different from the four models studied in both unitarised and 
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Figure 1: The e^e~ A*^/^" cross section. The data are shown to- 
gether with the fit from unitarised VMD2. Other models give similar 
descriptions. 



non-unitarised forms. 

The previous section illustrates that essentially all of the known formulations of VMD 
are able to provide a quite satisfactory description of the e~^e~ data. This was a conclusion 
previously reached in Refs. with models referred to in that work as Mi, M2 and M3. 
As the e'^e~ — > tt+tt" cross section is not sufficient to differentiate between the various 
realisations of VMD, comparison with other data and/or information is helpful. 

The e^e" tt+tt" cross section depends only on \F.,,{q'^)\, i.e., the magnitude of 
Ft^. We have at our disposal a number of models which give equally good fits to this 
cross-section. Hence we can now examine the ability of each of our versions of the VMD 
hypothesis to reproduce the I = I = 1 nn phase shift data, as we have arg [F^(g^)] = 6}, 
as noticed in Section p75| . 

As previously noted the resonant part of the form factor also contains in general some 
small isospin zero contributions {e.g., from the u and (p). Below 1 GeV, the dominant 
isospin 1 contribution to strong interactions is the p°. This was the conclusion previously 
reached when analysing vrvr data |38-40|. Moreover, it has been shown recently in 



r decay to pion pairs, that the infiuence of possible higher mass vector isovector mesons 
(p(1450) and p(1700)) is negligible below ~ 1 GeV. As r decays and e~^e~ annihilation 
are connected through CVC which has been shown in the conclusion of aleph 
gives support to neglecting these states in our invariant mass range. 

Thus, there is an apparent elasticity of the P-wave isovector vrvr scattering up to 
~ 900-f-lOOO MeV. Then , we can limit the relevant resonant contribution to the 1 = 1 = 1 
TTTT phase shift to the single p meson propagator. For VMD2, this means that the opposite 
of the p°-propagator phase coincides with the phase 5| of the I = I = 1 tttt partial wave 
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Figure 2: e~^e tt+tt cross section. The data are shown together with 
the fit curve obtained using VMD2 not unitarised. 

(see Eq. (0))- When deahng with all other models, there is a nonresonant direct ^tttt 
contribution to the pion form factor which also contributes. It is the concern of the 
present section to check to which extent, the phase of the form factors defined in Section 
H (together with the numerical values obtained from fit to e^e^ tt+tt^ and listed 
in Tables |1] and 0) gives an appropriate description of 6l, as it is determined and/or 
measured in nn scattering. 

TTTT scattering is an important physical process, where the basic concepts of S-matrix 
theory can be applied with some precision; these basic requirements include unitarity, 
analyticity and crossing symmetry. In the case of the vrvr system, crossing symmetry puts 



stringent constraints on the partial waves by mean of the Roy equations or dispersion 



relations. Methods based on these general principles are well known (see, e.g., Ref. 
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Figure 3: e^e tt+tt cross section. The data are shown together with 
the fitted curve obtained using unitarised VMD2. 



and references therein) and have been applied |13| in order to fit and/or reconstruct the 
1 = 1 = 1 TTTT phase shift 6l at low energy {i.e. from threshold up to about 1 GeV). These 
methods and most of their results are presently encompassed by ChPT |]5ll-p4|. 

Therefore it is of some interest to compare the phase shifts predicted from the versions 
of VMD obtained after fitting to the e'^e" vr+vr" data to the values of 6l as tabulated 
by Ref. [|l^] (see their Table 1). In this way, it is possible to check the consistency of the 
information deduced assuming each VMD formulation with the results of Ref. [|I^ which 
were derived under completely independent assumptions (namely, unitarity, analyticity 
and crossing symmetry). One should note that our model errors (those produced by the 
errors on the fitted parameters) have a small effect on the predicted phase; indeed, the 
accuracy on the predicted phase is 1.5% at 600 MeV and 0.6% at 800 MeV. From a 
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practical point of view, Ref. [|r^ does not quote errors. Therefore, we shall unfortunately 
not be able to express this comparison in terms of statistical quality factors. 
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Figure 4: tt+tt^ phase shift 6\. In both figures, the points are those from 
Ref. [|^. In (a) curves displayed correspond to non-unitarised models 
(parameter values taken from Table 0); full line is VMDl, dotted line is 
HLS and dashed line is VMD2. In (b) the curves displayed correspond 
to unitarised models (parameter values taken from Table the upper- 
most curve (dotted) is HLS, the lowest curve (full line) is VMDl; in 
between, the two curves corresponding to VMD2 and WCCWZ almost 
superimpose at this scale. 

In Fig. ^, we have plotted the curves predicted from our analysis and the phases 
from Ref. [0. The non-unitarised versions of VMDl and HLS can clearly be rejected 
as failing to support the identification of the b\ phase shift with the phase of ^^^(g^). 
Non-unitarised VMD2 looks in better agreement with expectation; however, recall that 
non-unitarised VMD2 is in poor agreement with e^e~ data (see Table above). These 
remarks are even clearer from Fig. ^, where the difference between the model phases and 
the data point phase of JTSf is plotted. Hence, we conclude that all non-unitarised forms 



of VMD are disfavored by the data. 

Instead, Fig. ^ exhibits remarkably good agreement of all unitarised versions of VMD 
with the data from Ref. |T3] over the full invariant mass range (from the two-pion thresh- 



old up to 1 GeV). Taking into account the simphcity of our expressions for ^^^(g^), this 
is already remarkable. The difference between the various predictions for the phase shift 
and the data of Ref. |13[ is presented in Fig. ||(b) for the unitarised models. Clearly, 



VMDl fails to describe Froggatt-Petersen phase over the whole invariant mass range, at 
a finer level of accuracy. Both WCCWZ and VMD2 have a systematic disagreement with 
the Froggatt-Petersen phase, (it is larger at low energy for WCCWZ than for VMD2). 
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Finally, the HLS model exhibits a quite remarkable agreement with the expected phase 
up to about 700 MeV, i. e. where the Froggatt-Petersen reconstruction is expected to be 
highly accurate; in this region the distance to expectation grows slowly, up to only 0.5 
degrees near 700 MeV. 

Given that the pole location for the p in the complex s-plane is an input into the 
phase shift analysis^ (set at a mass of 767 MeV and a width of 137 MeV and not allowed 
to vary), the disagreement of about 6 degrees near the p mass for the HLS model should 
not be considered presently as particularly significant and, correspondingly, the fact that 
the point where the phase goes through 7r/2 in Ref. (780 MeV) is closeQto the mass 
given by VMD2 (780.4 MeV) should also not be regarded as especially remarkable^- 

Finally, the relative difficulty for our unitarised models to match the phase above 
800 900 MeV (see Figs. ^3 and ||b) can also be partly attributed to higher vector 
mesons which definitely exist 0,^] and may indeed contribute at a small level in this 
mass region. However, such additional contributions should not be significant below the 
p peak. One should also notice that the phase shift of |T3] is not expected to be as 



accurate in this invariant mass region as it should be below the p mass peak. For all 
these reasons, the matching of phases in the region 800 -i- 900 MeV cannot be considered 
as constraining as the matching in the region 300 ^ 600 MeV. 

Another way to carry on the comparison of the FT^{q^) phases with the Froggatt- 
Petersen phase shift, is to superimpose the various model predictions for sm5\/p^. This 
function is connected with the isospin 1 P-wave scattering length a\ through: 

and has the characteristic of highly magnifying differences due to the phases in the low 
energy region (up to, say, 600 MeV). The curves corresponding to the various unitarised 
models are shown in Fig. ^ together with the points reconstructed using the phase as 
given in Ref. |T3|. Note that erratic behaviour of the values deduced from Ref. at the 



very beginning of the curve is entirely due to rounding errors only|^. The general success 
of the HLS model in matching the low energy tttt phase data is especially remarkable. In 
terms of angles, in the low energy region the distance to the data points does not exceed a 
few hundredths of a degree. Instead VMDl finds 0.02°, WCCWZ 0.03° and VMD2 0.06°, 
where the prediction is rather 0.1°. In Fig. |^ we present, for all unitarised models, the 
measured cross section for e^e" tt^tt^ and the previous fit function behaviours in the 
low energy region. They are equally as good with the possible exception of VMD2. From 

-"^^We thank CD. Froggatt and J.L. Petersen for this information. 
^"'^This is not an input in their method, but an output. 

^^In order for comparison at this level of refinement to be meaningful, the Froggatt-Petersen approach 
should first be carried out with a more accurate p pole location, those of Ref. |^ for instance. 

i^Ref. [|13| gives for the phases 0.1°, 0.4°, 0.7° at respectively 300, 320 and 340 MeV. In order to get 
a smooth behaviour consistent with the rest of the curve, it is enough to change these angle values to 
respectively 0.14°, 0.36°, 0.66°. This remark illustrates the magnification effect produced by the function 
in Eq. (pOD. 
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Figure 5: tt+tt phase shift S}. The plots display the difference between 
the phase shift of Ref. [0| and our phase shifts. In (a) curves displayed 



correspond to non-unitarised models (parameter values taken from Table 
full line is VMDl, dotted line is HLS and dashed line is VMD2. In (b) 
the curves displayed correspond to unitarised models (parameter values 
taken from Table |^); uppermost curve (dotted) is HLS, lowest curve (full 
line) is VMDl and middle curves are (dashed line) VMD2 and (dashed- 
dotted) WCCWZ . 



this, we conclude that the difference in ability to describe the 6\ phase shift is due to 
the existence and the magnitude of a direct ^TT^n~ couphng within models when fitting 
e'^e~ data. 

The success of the HLS model in reproducing the correct phase values and behaviour 
is remarkable. This can be partly attributed to the fact that its 77r7r coupling is connected 
with universality violation, while for VMDl, this term is fixed in absolute magnitude to 
the pion charge. Another reason is connected with the fact that in the HLS model, Qp^ 
is constant while it is proportional to in VMDl and WCCWZ. This actually produces 
a too large drop of the resonant contribution to the phase at low invariant mass which is 
not supported by the data, as it is illustrated in Fig. |. 

On the other hand, the WCCWZ model, while differing from all other models by 
having a conventional p mass and an additional dependent term, does not match the 
phase shift better than VMDl or VMD2. We may therefore conclude that the extra term 
in WCCWZ is not consistent with appropriate low-energy behaviour when WCCWZ is 
constrained to fit the broad energy range (0.28 to 1.0 GeV) rather than near threshold 
only. 
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Figure 6: Function sin 5 J /p^ deduced from -n^-n phase shifts; the func- 
tion is given in units of m~^. The dots correspond to the points of 
Ref. 0. Full line is VMDl, dashed line is VMD2, dotted line is HLS 
(the best description of expected phase) and dashed-dotted is WCCWZ. 



Figs. I, H and ^ demonstrate that, within the class of models considered, a good fit to 
a// experimental data can only be obtained by including an explicit direct 77r+7r^ coupling 
term in the model. For the HLS model, the magnitude of this direct coupling is quite 
well- fixed by data and is approximately — e/6 (in contrast to the common form of VMD, 
VMD2, for which such a term is absent). 

Indeed, Fig. |^, clearly shows that the four unitarised models fit practically as well 
the cross section e~^e^ tt+tt^ in the low energy region. Therefore, their relative ability 
to account here for the phase 6l does not originate from sharp differences in fitting 
e^e^ TT+TT" in the same region. The phase corresponding to the HLS model matches 
accurately 6l from threshold to about 700 MeV, as illustrated by Figs. § and |; 
this phase coincides with that of F-^i^q"^). Instead, VMDl having a 77r+7r~ coupling of 
magnitude fixed to e, fails to describe accurately this low energy region. VMD2 exhibits 
a systematic disagreement, which can be traced back to the lack of an explicit 7'7r7r 
coupling, more important at low energies; it should be remarked that VMD2 is, however, 
the closest (after HLS) to the data points; this is related to the fact that VMD2 actually 
corresponds to the HLS model with a = 2 while the data prefer a slightly different value 
(a ~ 2.4). 

The behaviour of WCCWZ with respect to HLS clearly shows that data considered 
from threshold to about the mass are more consistent with the existence of a constant 
nonresonant coupling. Another interesting result is that, from Table 0, it is clear that we 
can perform the fits with the HLS model fixing A = 1 without any significant change in 
the fit quality or in the parameter values. The fact that A ~ 1 is preferred by the data 
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at the fit level is interesting in a number of regards; indeed, A = 1 is one of the most 
common ways to parametrise the p width; it is also the value obtained for the width for an 
interaction term in the Lagrangian of the form p6| p^(7r+9^7r^ — 7r^9^7r+) with constant 



dressed coupling. We discuss in Section whether a departure of A from 1, even if not 
significant at the fit level (see Table ||), should be considered. In closing this section we 
note that no non-unitarised model appropriately reproduces the nn phase shift |]TB[ , and 
that, among all unitarised models we have examined, only the HLS model of Bando et 
al. agrees well with the tttt phase shift predictions [O . 



6 ChPT and constraints in the near-threshold region 

In the preceding sections, we have analysed various realisations of VMD using e^e~ data 
and a tttt phase shift analysis, over the full invariant mass range where data and/or 
information is available. It is also of interest to examine in detail the consequences of the 
various VMD models for the magnitude of -F^r in the near threshold region. 

From a general perspective (such as that of ChPT, for example), the VMD La- 
grangians are rather simple and should not be expected a priori to be able to provide an 
accurate representation of the full physics of QCD over the whole of the low-energy regime 
below 1 GeV^. In particular, it is known that the effects of p exchange first show up in 
the usual Gasser and Leutwyler effective chiral Lagrangian, CeE, (relevant to processes 
involving only external vr legs as contributions to the fourth order low-energy 



constants (LEC's), L^, appearing in C^s This means that, very near threshold, other 



contributions, not generated by p exchange, are also, in principle, present in vrvr scatter- 
ing and Ft,. The models VMDl, HLS and WCCWZ allow for such a possibility, though 
in a form much simpler than would necessarily be expected on general principles. The 
fact that we allow some phenomenological freedom in the treatment of the s-dependence 
of the p width, however, means that some of the more complicated s-dependence of the 
non-resonant contribution may be able to be reproduced by a slight adjustment of the 
parameter A. Such a readjustment, of course, has no effect at the p peak. 

The physical spectral function near threshold in the vector isovector channel, which 
is measured experimentally in e^e~ tt^ti^ and in r decays [Q, is very well described 
near threshold by ChPT |]55|, |5^, ^ . At next-to-leading (two-loop) order in the chiral 
expansion, the ChPT expression for F-^lq"^) however, contains two new sixth order 



LEC's not constrained by other experimental data, so an improved prediction is not 
possible at present. However, one can use the form of the two-loop prediction |]59|, 



F.(g2) = 1 + UrXq' + c^q' + /^(gVO + 0{q'), (31) 
o 

(where fy is a known function, whose explicit form may be found in Ref. ||6^) to analyse 
the near-threshold data and obtain an optimal model-independent fit for the parameters 
(r^)y and c-,,. The results are [0] 

(r^)^ = 0.431 ± 0.020 ± 0.016 fm^ 
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= 3.2±0.5±0.9 GeV"^ 



(32) 



where the first error is statistical and the second error theoreticaL Note that, although 
the data analysed is that of the NA7 collaboration 
a factor of 2 larger than quoted in Ref. 



61[] , the statistical error is roughly 
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According to the authors of Ref. ||6^ this 
larger error reflects the presence of the additional fit parameter, c,^. We will accept the 
results of the analysis of Ref. since the form of the fit function used, being obtained 
from ChPT, is the most general one possible, compatible with QCD, to this order in 
the chiral expansion. As such, the analysis is as model-independent as possible, to this 
order. Using these results we obtain an experimental value for the threshold value of 
Fjr. In quoting this result below we have added all errors in quadrature. Note that the 
next-to-leading order contributions appearing in this experimental fit are in agreement 
with those appearing in the 1-loop ChPT expression for FT^{q^) [^. This is of interest 
since, in the case of the one-loop expression, the relevant LEG, Lg, is also constrained 
by other experimental data. Although it is conventional to fix Lg using (r^)y as input, 
fixing it independently through a quantity like the combination r^//iy appearing in the 
amplitude for the process 7r+ e'^u^e^e' leads to compatible values. As a result, the 
1-loop expression for may be thought of as a prediction. 

We now compare the results of the various unitarised models, as fitted to e~^e~ cross- 
section data, for F7r(4m^), with the experimental value, obtained as described above[3 



' [i^.(4m2 


)]vMDl 


= 1.174 ±0.001 




)]wccwz 


= 1.149 ±0.010 


< [FA^ml 


)]hls 


= 1.176 ±0.001 


[FA4^ml 


)]vMD2 


= 1.393 ±0.001 


[F^{4:ml 


)]Expt 


= 1.17±0.01 



(33) 



We see that each of these unitarised VMD models is close to the ChPT result with 
the exception of VMD2. 

In order to quantify our comparison of the low energy phase behaviour we have com- 
puted the scattering lengths a\ following from the various unitarised models. For this 
purpose, Eq. ([30|) has been used with the corresponding phase of each F.„, following our 



remarks in Section 2.5; the results are 



' ai(VMDl) 
a\(VMD2) 
ai(HLS) 
ai(WCCWZ) 



0.006 ±0.007 
0.023 ±0.002 
0.043 ±0.003 
0.013 ±0.002 



(34) 



^"'There was a small change in the "experimental" number above from the fy term in the general 
expansion. 
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in units of m^^. These values can be compared with experimental results from Ke^ 
data |63| using a Roy equation fit (aj = 0.038 ± 0.002). The result of Ref. 



[a 



0.10) disagrees with these results and that treatment violates basic assumptions of 



ChPT as illustrated in Ref. [Q. The Current Algebra prediction |]65| is a\ = 0.030 and 
the ChPT result |5| predicts a\ = 0.037 ± 0.001 at the two loop order (at 0(j9^)). A 
recent preliminary ChPT calculation predicts that the expected value for a\ should 
increase at order 0{p^) and give a value in the range 0.038 0.040. The range given for 
this prediction reflects the existing uncertainties in the low energy constants (LEC's). 



The question with the above estimates of the scattering length a\ is twofold : 1) to 
what extent is the interpolation of each model into the low energy regime and subsequent 
extraction of a\ meaningful and reliable? 2) is it possible to improve each estimate? 

Assuming that the interpolation to low energies is meaningful, which seems legitimate 
in view of the preceding section (at least for the HLS model), the answer to the "reliabil- 
ity" issue is connected with the fit quality for each case, which is well summarised by the 
line x^/dof(7r7r) in Table ^. As noted earlier, the best fit quality reached with the existing 
e~^e~ —>■ Ti^-K" data corresponds to x^/dof = 61/77. This value was already achieved 
in Refs. with their model Mi; we get it here with the WCCWZ model which is a six 
parameter fit while Mi was a five parameter fit (as our present models VMDl, VMD2, 
HLS). From a statistical point of view, a x^/dof = 61/76 is unlikely to be reliably im- 
proved and should be considered close to perfection and hence, aj (WCCWZ) cannot be 
significantly improved. Taking into account the systematic disagreement with the phase 
of Ref. |TB[ revealed by Fig. ^ we cannot consider this estimate as reliable. 

Unitarised VMD2 gives a very good (even if not perfect) fit to the pion form factor 
(x^/dof = 81/77); however, from Fig. |^, it appears to have some difficulty accommodating 
all data points as well as the three other VMD models in the low energy region; this can 
be attributed to the lack of a non-resonant term which is present in the other models. 

Finally, VMDl and the HLS model provide both a near optimum fit (x^/dof = 65/77) 
to the pion form factor, and so cannot be dramatically improved. We thus believe that, 
numerically, the extrapolation is reliable and cannot be significantly improved. We shall 
nevertheless examine the consequences of fixing A = 1 in the next section. 

To summarise this section, we can say that among the unitarised versions of the 
VMD assumption, each of VMDl, WCCWZ and HLS provide a good extrapolation for 
the magnitude of all the way down to threshold, while only the HLS model provides 
additionally a good extrapolation of the phase 5\. From the discussion above, we can 
give our extracted value for the scattering length, a\ = 0.043 ± 0.003, which compares 
well with the previous measurement [^,^ (their separation is 1.4o"). 

As a consequence of these observations we conclude that, of the models considered, 
only the unitarised HLS model, with a ~ 2.4, provides a successful representation of 
in both magnitude and phase which is valid in the whole region from threshold to the (p 
mass. 
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7 Discussion 



The previous analysis of the e~^e data and its comparison with the phase shift 
51 of Froggatt-Petersen [0], lead us to conclude that the VMDl, VMD2, WCCWZ and 
HLS unitarised models provide a good description of a very large set of experimental 
data, as illustrated by Table |^ and by Fig. ^. However, a finer check on the low energy 
behaviour (see Fig. ^ highly favors the HLS model over all other ways to formulate the 
VMD assumption. This conclusion is enforced if one includes the known information in 
the near-threshold region provided by experiments and ChPT (values for Fj^{4'm'^) and 



This agreement is reached by allowing the HLS a parameter to slighly depart from 
its universality value: a ~ 2.4 instead of a = 2. Relying on this model, we find that 
universality is violated at a level of 20%. As a direct consequence of this, we observe a 
non-resonant direct coupling 771"'' vr", with strength —0.182 ± 0.008. The corresponding 
term in VMDl is 1 (due to the pion electric charge), while it is for the widely accepted 
VMD2. Universality violation does not prevent the HLS model fulfilling the constraint 
F^(0) = 1 in a quite natural way. 

Even if the influence of this direct term is the largest near threshold, it contributes a 
significant improvement up to 1 GeV, because of interferences. This becomes even clearer 
if one notes that VMD2 is qualitatively equivalent to the HLS model by simply removing 
the constant term^ in its expression for (see section 0). Therefore, we have the first 
suggestive evidence for a non-zero point-like 77r^7r^ coupling. 

The origin of the shortcomings of the VMDl, VMD2 and WCCWZ models are read- 
ily understood. For VMD2, universality violation (required by the data) produces an Ft, 
which, extrapolated to = 0, has F.jr{0) — 1 = 6 > 0. Not surprisingly, near threshold 
then, the fitted VMD2 model exceeds data by almost exactly S. For VMDl and WCCWZ, 
the problem is the near-threshold phase. In a conventional effective field theory formula- 
tion ||18||, however, the direct pion-photon coupling term in (unavoidably present) has 
a non-zero phase via vrvr rescattering. The truncated version represented by VMDl and 
WCCWZ does not allow for this possibility and the models turn out to be insufficiently 
flexible to fully compensate for this deficiency by an adjustment of the parameters of the 
p contribution. 

An interesting issue is also the observed inequivalence of VMDl and VMD2. One 
may argue that this could be a shortcoming due to having neglected mass dependence in 
the p — 7 coupling produced by loop effects. Following Ref. p6|, we have checked that 



this is not the case, mainly because universality violation is clearly preferred by the data. 

Another interesting issue is the value found for A (1.06 ± 0.08) in the unitarised 
HLS model. The departure from 1 is not statistically significant as far as fits of e'^e~ 
data alone are concerned and hence the question of whether one can set A = 1 should be 
considered. Because of unitarisation (see subsection |2]^) , there is no reason why one could 
not observe departures from 1, which could effectively account for neglected contributions 
{t— and u— channel exchanges, non-resonant strong interaction among pions, higher mass 

^^The role played by a in gp^ in the HLS model is then transferred to e in VMD2. 
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resonances, etc.). The fit value for A sliows tliat tliese neglected effects are small enough 
not to spoil the analytical shape generally expected for the mass-distribution as is 
obvious from Table 0. However, if all consequences of fixing A = 1 were acceptable, this 
will allow us to provide a 4 parameter fit to the pion form factor (the smallest possible 
set): 3 of these parameters describe the p mass and couplings to ti^ti~ and e"'"e~, and the 
Orsay phase which is tightly connected with the omega isospin violating contribution. 

We have done the fit with the HLS model, fixing A = 1 in order to get the full set 
of parameters needed to check accurately the near-threshold results of the preceding 
section. We have also left free the parameter A which governs the magnitude of Br(co' — > 
TT+TT"), as it is mainly influenced by e^e~ — tt+tt" data. As expected, the fit quality is 
practically unchanged (x^/dof = 66/77). We have obtained for F7r(4m^) = 1.177±0.001, 
in good agreement with the ChPT expectation and the value reconstructed from other 
experimental data [§^,0 (1.17 ± 0.01). We have also obtained a\ = 0.041 ± 0.003 using 
the full phase of Fj^{q^), which is in better agreement with the two-loop calculation of 
(0.037± 0.001). Moreover, our result is in fairly good agreement with [p^] which predicts 



ai = 0.040 ± 0.003, relying on Roy equation techniques. This is consistent with the vrvr 
P-wave being nearly purely resonating with only a small background from threshold to 
900 MeV. It therefore appears legitimate to fix A = 1 in our fitting procedure with the 
unitarised HLS model. 



We give as final results those obtained in this last kind of fit with the HLS model |17 



mp(MeV) 


= 775.1 ±0.7 


r,^,+,-(MeV) 


= 147.9 ±1.5 


Fp^e+e- (keV) 


= 6.3 ±0.1 


Bt{uJ TT+TT^) 


= (2.3 ±0.4)% 




= 1.177 ±0.001 


<P 


= 0.041 ±0.003 
= 104.7° ±4.1° 


a(HLS) 


= 2.37 ±0.02 


A 


= 1 (fixed) 



(35) 



Concerning F^(4m^) the value above is in good agreement with ChPT predictions 
and the previously quoted experimental value. Our value for a\ compares well with the 
determination given in which relies on K^.^ data and the use of Roy equations ||49|| , 
and with the ChPT calculation of |53] and the Roy equation result [^|. Taking into 
account the uncertainties on the low energy constants (LEC's), ChPT at order appears 



consistent with a\ = 0.040, in good agreement with our extracted value. 

A final remark concerning the HLS model is of relevance. As said above, it predicts 
that the strength of direct •yini coupling (let us call it provisionally c) is tighly connected 
with universality violation (and then to a small violation of the KSFR relation) by 
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c = 1 - ^^^^ (36) 
ml 

This relation is found to give c = —0.182 ± 0.008 when the constraint c = 1 — a/2 
is forced. As a matter of check, we have rerun our fit procedure by decouphng c from 
Qp-y and QpTTTT- This corresponds exactly to the model M3 of [§, however leaving free 
gp-y as it should. We have, of course, obtained a very good fit (x^/dof = 66/76), fixing 
also A = 1, which provides c = —0.171 ± 0.018, in fairly good agreement with the HLS 
model result]^. In other words, the HLS model is able to provide a meaning to the 
direct coupling c in terms of only the parameters (see Eq.(PBD just above). This (non 
obvious) correlation was, of course, completely missed in and can be considered as a 
remarkable success of the HLS model. Finally, the failure of VMDl and WCCWZ tends 
also to support the conclusion that Qp^ is more consistent with a constant value, than 
carrying a dependence. Indeed, this dependence produces a too strong suppression 
of the p contribution when going down to threshold. 

The question is now how to compare the other p° parameters given above to the 



corresponding existing measurements |41|. The main problem with an object as broad as 
the p^ meson is that its shape and, correspondingly, its observed parameters are highly 
influenced by phase space effects and by the production mechanism (creation amplitude). 



This explains the wide spectrum of values reported in the Review of Particle Physics ||41[ 



Generally, phase space effects are perfectly known and departures from expectations 



have a physical meaning, as in the decay rj' 7r+7r~7 (see Ref. P, 37| for instance). 
In hadronic and in photoproduction experiments, one has to rely on expressions for the 
production amplitudes tcN pN and — > pN, which are actually guesses to a large 
extent; Ref. [^] illustrates the dependence of the p mass on several (and all reasonable) 
guesses for the vrA^ — >■ pN amplitude. In the case of e^e~ annihilations, we also work 
under assumptions which can be discussed (see Ref. [|,|| for instance) mainly about the 
transition 7 — p° and the possible existence of a constant 77r+7r^ coupling. One of the 
original motivations of the present work was indeed to rely on both e~^e~ vr+vr" and 
e~^e~ p^p~ data in order to get rid of a possible 77r+7r~ coupling. We have shown 
that the single existing data set with pp final state is not precise enough in the p^ region 
in order to achieve this program. However, it happens that the existing information on 
6l, obtained in a completely independent way, allows us to single out the influence of 
the 'J p^ transition and those of the 77r+7r~ coupling. The quality of the comparison 
between the phases of the HLS model on the one hand and the data points of Ref. [|T^ 
on the other hand shows that these two sources of errors are well controlled. 

The fit values obtained for our free parameters and given in Eq. (^3p are calibrated by 
a fit on e^e' tt+tt^ data, i.e. they optimise only |F^(g^)|. The expression deduced for 
arg[F7r(g^)] automatically matches all known information on the phase Si without any 
further tuning, as illustrated in Sections || and |[ Therefore, we conclude that our results 
in Eq. (RSI) are little affected by systematic errors due to modelling. The single remaining 



In with gp-y fixed at its PDG value and A left free, the result was c — — 0.26± 0.04; with this 
respect, note a misprint in the caption of Table 1 here: A = co + C2m^ should be read A = —{co + C2m^). 
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freedom in defining the p° parameters is the definition used for mass and width. What 
we have quoted corresponds to the most usual Breit-Wigner definition. 

Moreover, ChPT predictions |]53|,|56[| near the two-pion threshold allow us to check 
even more accurately the quality of the information deduced from e~^e~ data. We can 
therefore conclude that the HLS model carries information reliable enough to lessen to a 
large extent the dependence of the parameters on the production mechanism. 

An improvement over the results given in Eq. (|35D would require comparison with 
other independent data and an updating of the Froggatt-Petersen spectrum |TB|, using 



more recent data and a more accurate p pole position |^ 

and e 



Moreover, the accuracy of 
is limited by large systematic 



the existing data on e~^e~ — > 7r"'"7r ana e'e — /i ' /i 
uncertainties. One can expect a dramatic improvement from the new experiment with 
the cmd2 detector, now in progress at Novosibirsk, aiming to reach a statistical accuracy 
of 3% per point and an overall systematic uncertainty better than 1% p8| . 

Finally, our results suggest a HLS-like model might also prove useful in parametris- 
ing hadronic spectral functions measured in r decay where, to this point, VMD2-like 



parametrisations have typically been used [43 



8 Summary and Conclusions 

We have studied a variety of vector meson dominance models in both non-unitarised and 
unitarised forms. They depend on a few parameters (mass of the p meson, its coupling 
constants to vrvr and e~^e~ , the shape parameter A and the Orsay phase needed in order to 
describe the p—u mixing). We have fitted these to both e~^e~ — > tt+tt" and e~^e~ — > p~^p~ ■ 
In order to study the behaviour of each solution, we have studied how they match the vrvr 
phase shift obtained under general model independent assumptions from threshold up to 
1 GeV. We have also examined the value they provide for threshold parameters [F^^ at 
threshold and the scattering length a\), which can be estimated accurately from ChPT. 

This represents the largest set of independent data and cross-checks done so far. It 
happens that, of the models considered, only the unitarised HLS model is able to account 
for all examined effects. We also find that the standard value A = 1, corresponding to a 
point-like pnir coupling, is well accepted by the data for p parametrisation. 

Unlike the standard formulation of VMD, fits with this model return a significant 
non-resonant contribution to the electromagnetic pion form factor. This was found to 
have a value ~ — e/6. This term is governed by a small universality violation which 
changes the HLS parameter a from 2 to 2.4. All other models considered, even if they 
are able to describe e~^e~ annihilations quite well, are unable to account satisfactorily for 
the other available information. It should be mentioned that, within the class of models 
considered, our results tend to favor a constant gp-y over a dependent one. 

We give the values for the p° mass and for its partial widths to 7r~^7i~ and e~^e~, 
obtained using the HLS model with A = 1. This corresponds to describing the resonance 
mass spectrum through the usual Breit-Wigner expression (A = 1). We also obtain an 
estimate of Br(ti; tt+tt^). 
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As such, we conclude that, of the models considered, the HLS model with a ~ 2.4 is 
the most favoured version for implementing the VMD ansatz. Thus, it is interesting to 
consider whether the success of its predictions for the magnitude and phase of F^(q'^), 
could be obtained in a way which does not need the assumption that the p is a dynamical 
gauge boson of a hidden local gauge symmetry. Hence, looking for other models able to 
describe, as successfully as the HLS model, the same large set of data is useful in order to 
know whether the conceptual motivation for this model should be interpreted as having 
any underlying significance. 
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Figure 7: e^e 7i~^n cross section in nb with fits using various uni- 
tariscd VMD models. Data point and curves are given from tliresliold to 
500 MeV. 
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